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ABSTRACT 

SN 1997D in NGC 1536 is possibly the least luminous and energetic Type 
II supernova discovered to date. The entire light curve is subluminous, never 
reaching My = —14.65. The radioactive tail follows the 56 Co decay slope. In 
the case of nearly complete trapping of the 7-rays, the 56 Ni mass derived from 
the tail brightness is extremely small, ~ 0.002 M . At discovery the spectra 
showed a red continuum and line velocities of the order of 1000 km s -1 . The 
luminosity and the photospheric expansion velocity suggest that the explosion 
occurred about 50 days before discovery, and that a plateau probably followed. 
Model light curves and spectra of the explosion of a 26 M star successfully 
fit the observations. Low mass models are inconsistent with the observations. 
The radius of the progenitor, constrained by the prediscovery upper limits, is 
Ro ^ 300 R Q . A low explosion energy of ~ 4 x 10 50 ergs is then required in 
the modeling. The strong Ba II lines in the photospheric spectra are reproduced 
with a solar abundance and low T e g. A scenario in which the low 56 Ni mass 
observed in SN 1997D is due to fall-back of material onto the collapsed remnant 
of the explosion of a 25-40 M & star appears to be favored over the case of the 
explosion of an 8-10 M & star with low 56 Ni production. 

Subject headings: nucleosynthesis - stars: evolution - supernovae: general - 
supernovae: individual (SN 1997D) 
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1. Introduction 

SN 1997D was serendipitously discovered on Jan. 14.15 U.T. ( pc Mcllo fc Bcnctti| 



1997|) in NGC 1536, a morphologically disturbed spiral galaxy belonging to a high density 
group QMaia et al. 1994Q . Although NGC 1536 {v hc i io = 1296 km s"\ RC3) is one of the 
galaxies patrolled visually by Rev. Evans for his SN search (private communication), he 
missed the SN because its brightness apparently never exceeded his detection limit. 

The first spectrum showed the main features of type II SNe, but also revealed 
peculiarities which make this object unique. Particularly noteworthy were the extremely 
slow expansion velocities, the red color and the strong Ba II lines ( |Turatto et al. 1998|) . A 
campaign of photometric and spectroscopic observations was therefore promptly started at 
ESO and CTIO. The complete data set will be presented elsewhere. 



2. Properties and evolution 

The SN was intrinsically faint compared to other SNe II. If we assume a distance 
modulus of /j, = 30.64 mag for NGC 1536 (H = 75 km s" 1 Mpc" 1 ; |Tully 1988Q , and 



Ab = 0.00 (RC3), a set of closely spaced non-detections by Evans sets an upper limit to 
the SN luminosity of Mv(max) £ —14.65. This is about 2 magnitudes fainter than regular 
SNe II at maximum ( IPatat et al. 1994] ). 



In Fig. [TJwe compare the absolute V light curve of SN 1997D with those of other SNe II: 
the double-peaked SN 1987A, the prototypical SNIIb 1993J and the linear SN 1979C. Even 
with the uncertainties on the epoch of maximum and in the distance modulus, SN 1997D 
was definitely fainter than all of these objects. The decline rates in the V, R and I bands 
after J.D. 2450520 are close to the 56 Co decay rate. 

The color of SN 1997D at discovery was unusually red. No sign of interstellar absorption 
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lines, which are expected in case of strong extinction, is present in the spectra. SN 1997D 
reddened rapidly over the first few weeks after discovery, reaching a maximum (B— V)= 2.5 
in mid-February, and then turned gradually bluer. The analogy with the color curves of 
other SNe II, which reach the red-most point between 40 and 80 days after peak luminosity, 
indicates that maximum light should have occurred around J.D. 2450430 (±20d). 

A bolometric light curve is needed for a proper analysis of the energetics. Unfortunately, 
most of our observations only cover the optical domain. Only at one phase are nearly 
simultaneous optical (Jan. 31) and IR (Jan. 27) spectra available. The overall luminosity 
between 0.35 and 2.3 \i in the combined spectrum is log(Lboi) = 40.83 erg s _1 . At that epoch 
the color is still sufficiently blue that Schmidt's (1998) empirical bolometric corrections 
to V can be used. One then obtains log(Lboi) = 40.96 erg s _1 , in good agreement with 
the spectrophotometry. At later phases the color of SN 1997D was so red that Schmidt's 
correction cannot be applied. However, the available spectra show that at all epochs most 
of the energy was radiated in the V and R bands and that the color evolution is mostly 
due to variations in the B band, where the flux is much smaller. Therefore, during the 
photospheric phase we assumed that the bolometric light curve can be described reasonably 
well by applying a constant bolometric correction to the V magnitudes. The bolometric 
light curve thus derived is shown in the insert of Fig. [TJ. The last point has been obtained 
adopting for SN 1997D the same L bol /L B VRi ratio as for SN 1987A in the nebular epoch 
( [Schmidt 1998]) . 

Although in SN 1987A the 7-rays from the 56 Co decay were not completely deposited 
QPinto et al. 1992]) , the late time decline of SN 1987A was close to the decay rate of 56 Co. 
The mass of 56 Ni ejected by SN 1987A was ~ 0.075 M Q QBouchet et al. 19914 |Bouchet et 



al. 1991 b| , [Bouchet fc Danziger 1993| , |Danziger et al. 1991| , [Nomoto et al. 1994| ). Assuming 



for SN 1997D the same 7-ray deposition as in SN 1987A, the ratio of the luminosities 
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on the tail of the light curves can be used to determine the mass of 56 Ni synthesized in 
the explosion of SN 1997D. We then obtain a ratio M 5 6 Ni (87A)/M 5 6 Ni (97D) ~ 40. Taking 
into account the uncertainties in the determination of the epoch of maximum and of the 
bolometric luminosity, the mass of 56 Ni ejected by SN 1997D is very small, between 0.001 
and 0.004 M . 

A comparison with SN 1987A (Fig. H) shows general similarities with respect to 
P-Cygni lines of H I, Ca II, Na I and Ba II, and more subtle differences. These latter 
are emphasized by the strengths of the Ba II absorption lines and possibly those of other 
s-process elements. SN 1987A showed strong lines of Ba II, the analysis of which ( |William"s 



1987 




Lucy I98*i, 


Hohich 1988 


J 



of barium. SN 1997D shows even stronger lines of Ba II, with the identification of an array 
of multiplets, but, as we shall discuss later, an overabundance of barium is not required 
because a lower temperature prevails in the envelope. 

Strong lines of Sr II, also an s-process element, were also observed in the IR spectrum 
taken on Jan. 27 ( piocchiatti et al. 1998] ) . In this spectrum the He I 10830A line, which is 



the dominant feature in the IR spectra of other SNe II, is not clearly detected, although it 
might affect the blue wing of Sr II 10914A. 

Finally, one of the most striking peculiarities of SN 1997D is the low expansion velocity 
of the ejecta, approximately 1200 km s _1 , as measured from the line absorption minima. 
Only the strong Ba II lines show somewhat higher velocities, but it is well known that 
strong lines may form in outer parts of the envelope and that stratification effects may also 
be present. As we will discuss later (Sect. 3.2), the small velocity constrains the epoch of 
the first spectrum to be ~ 50 days. 
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3. Model Light Curve and Spectra 

The light curve of SN 1997D is unique because the luminosity is low, both at 'maximum 
light' and on the radioactive tail. These properties, along with the narrowness of the lines 
and the redness of the continuum, can be used to deduce the physical characteristics of the 
progenitor and of the explosion. 

The progenitor model used here is that of a 26 M© star with an 8 M© He core at the 
beginning of collapse ( [Nomoto fc Hashimoto 1988|) and an 18 M© H-rich envelope. The 
stellar radius is i? = 300 R®, which is 3 - 5 times smaller than typical red supergiants. 
The explosion was simulated using a 1-D Lagrangian hydrodynamical code ([Young et al. 
1995[ ) with an energy of E = 4 x 10 50 ergs and the ejection of 0.002 M© of 56 Ni. In order 
to eject a small amount of 56 Ni, most of the material of the innermost region where Si is 
burned to Ni has to fall back on the degenerate neutron core where it is photodisintegrated 
before emitting 7— rays in forming the compact remnant. This is used to place the mass 
cut, resulting in a collapsing core of 1.8 M , and ejecta mass of 24.2 M©. 

3.1. Light Curves 

Our model light curve (Fig. |^) successfully reproduces the main characteristics of the 
observations. The model parameters are constrained as follows. The early light curve is 
powered by shock heating, i.e., the gradual release of deposited shock energy. After the 
peak, the model light curve enters a slight plateau phase lasting about 40 days (Fig. |3|); 
during this period, the recombination wave moves inwards in mass through the H envelope 
down to the He core. The peak luminosity {My £ -14.65) depends on Rq, M, and E. In 
order to achieve a consistent model, the progenitor radius cannot exceed 300 R & . For larger 
radii, the SN would have been too bright at shock breakout. 
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After the plateau the light curve drops over 2 magnitudes to the radioactive tail. 
The luminosity and shape of the tail depend on the 56 Ni mass and on E/M. The 56 Ni 
mass determines the absolute magnitude and E/M determines the decline rate. The 
tail of SN 1997D follows the 56 Co decay rate as in SN 1987A. For the SN 1997D model, 
E/M = 1.7 x 10 49 ergs M©" 1 . This is ~ 10 times smaller than in SN 1987A, and gives the 
low velocity and the complete trapping of the 7-rays. The low luminosity tail is reproduced 
with the low 56 Ni mass of 0.002 M . 

We have selected the current model through an elimination process involving 7 models, 
which we will discuss in a subsequent paper. For models with a small ejected mass (M 
< 20 M Q ) the plateau is too short, giving a time from explosion to discovery of < 50 
days, inconsistent with the spectra. The tail of SN 1997D falls at the 56 Co decay rate, 
which suggests two possibilities: 1) the ejecta is optically thick to the 7-rays and the tail 
reflects the total 56 Ni mass ejected; or 2) the ejecta is optically thin to the 7-rays and the 
tail is powered by annihilation of the trapped positrons. Since the positrons carry only 
3.7% of the decay energy, this latter case would imply a normal 56 Ni production. This can 
only be true if the opacity to 7-rays is less than one, i.e. if the total mass of the ejecta is 
M e j < 4n ^ 2 ■ A simple estimate shows that for this to be the case the mass of the ejecta 
would have to be less than the mass of ejected 56 Ni alone. Using for the tail the values: 
7-ray opacity, k 7 = 0.03 cm 2 /g, velocity, v = 700 km/s and time, t = 70 days, we get 
M e j < -^M Q . Thus the tail cannot be powered by positron annihilation and we therefore 
explore models in which the ejecta is optically thick to 7-rays. Models with larger energies 
give higher material velocities and broader spectral lines, which is again inconsistent with 
the observations. For models with larger E/M the light curve is too bright and the plateau 
is too short. Models with smaller progenitor radii have plateaus that are too short. To 
attain the same final kinetic energy, the shock, in a more compact star, expands the ejecta, 
and thus cools it, much more rapidly than for a star with a larger radius. The result is that 
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the photosphere moves inward faster, and the plateau is shorter. 

3.2. Spectra 

We used the explosion model described above to compute synthetic spectra for two 
epochs with a Monte Carlo spectrum synthesis code, which is based on that described in 
Mazzali Sz Lucy (1993) but has been improved to include photon cascades (Mazzali & Lucy 
1998, in preparation). 

On Jan. 17 the red continuum indicates a photospheric temperature T e g ~ 6400 K. 
The other input parameters are the expansion velocity (t> p h = 970 km s _1 ) and the epoch of 
explosion. Given the luminosity as deduced in Sect. 0, the required temperature can only be 
achieved with a rather large radius of about 6000 R G . Since i? p h = v p ht, this implies t ~ 50 
days for the first spectrum. Therefore, we deduce that the peak was missed by about 50 
days, that a plateau followed, and that the SN was first observed on the decline to the tail. 

The model spectrum (Fig. |j) matches the observed one well with the exception of some 
features in the red, where the model continuum is also somewhat too high, and of the B.j3 
absorption. Most of the narrow absorption lines are well reproduced. In particular, Ba II 
6497A is stronger than Ha. Since solar abundances were used, this is just a temperature 
effect, and no overabundance of Ba is required. The temperature is low enough for the 
Ba III/Ba II ratio to be small, which results in strong Ba II lines, while the same low 
temperature makes Ha weak. The narrow lines are clearly the result of a low explosion 
energy. The fact that the lines are so narrow allows a reasonably accurate determination 
of the mass above the photosphere. Values ranging between 10 and 2OM give acceptable 
spectra. The model adopted here has 14.5M Q above 970 km s _1 . Models with a smaller 
mass give lines too shallow if w p h = 970 km s -1 . 
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Synthetic spectra show that the requirement that t ~ 50 days on 17 Jan. is a rather 
strong one. If the explosion took place closer to discovery, i.e. t < 50 days, the observed 
expansion velocity implies a small radius. Together with the derived luminosity, this results 
in a higher T e g than required. In this case, the Ba III/Ba II ratio increases and the Ba II 
line intensities decrease. If on the other hand the explosion occurred earlier, so that t > 50 
days, the radius is too large and the temperature becomes too low. When the temperature 
drops below 5000K, lines of Ti I and Fe I absorb most of the radiation below 4500A and the 
synthetic spectrum changes dramatically. 

We also computed a synthetic spectrum of our selected explosion model evolved to 
match the 6 Feb. data. The fitting was satisfactory using an epoch consistent with t = 50 
days on Jan. 17. 

4. Discussion 

SN 1997D is among the least luminous Type II supernovae known to date. The low 
luminosity persists also in the tail, where the decline rate is close to that of a 56 Co-powered 
tail. SN 1997D is also characterized by a low temperature and a very low expansion velocity. 
These features are well reproduced by an explosion model with a relatively small progenitor 
radius (300 R & ), massive ejecta (~ 24 M ), low explosion energy (4 x 10 50 ergs), and small 
ejected 56 Ni mass (~ 0.002 M ). 

A similarly small amount of 56 Ni has recently been proposed for another SN II, 
SN 1994W ( |Sollerman et al. 1998|) , but this object was rather luminous and its tail declined 
faster than the 56 Co slope. 

The neutrino-heating mechanism of massive star explosions is not well understood, so 
the explosion energy and the amount of 56 Ni ejected as a function of progenitor mass are 
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difficult to predict from the hydrodynamical models. According to pre-supernova model 
calculations, stars more massive than ~ 25 M form an Fe core in a greater gravitational 
potential because of the significantly smaller C/O ratio and the consequently weak carbon 
shell burning ([Nomoto et al. 1993| ; [Woosley fe Weaver 1995|) . Therefore, if the efficiency of 
neutrino heating does not change with exploding mass, more massive stars tend to produce 
lower explosion energies, and they suffer from more fall-back of material onto the collapsed 
remnant. The mass cut is then placed further out. However, the progenitor must be less 
massive than the Wolf-Rayet progenitors (i.e., ^,30 - 40 M & ) because of the presence of an 
envelope with a mass of at least 15 M above the photosphere at day 50. 

Stars in the range ~ 10 - 25 M are likely to explode with typical explosion energies 
of ~ 1 x 10 51 ergs, and to produce 0.07 - 0.15 M Q of 56 Ni, as indicated by the brightness 
and light curve shapes of SN 1993J, SN 1987A and of Type Ibc supernovae (|Nomoto et al. 



|1995| ; [Young et al. 1995|) . Stars of 8 - 10 M Q are predicted to produce little 56 Ni and other 
heavy elements (|Nomoto 1984j ; [Mayle fe Wilson 1988|) . However, these stars are at the top 



of the AGB in the pre-supernova stages (Hashimoto et al. 1993j ), which is inconsistent with 



the small radius of SN 1997D. Also, synthetic spectra obtained from the low mass models 
have shallow absorption lines compared to the observations. 

In conclusion, the progenitor of SN 1997D was probably as massive as 25 - 40 M & . 
Our 26 M Q model is in this mass range, thus being consistent with the low E/M and the 
small 56 Ni mass. The presupernova radius depends on many parameters, and the reason 
for its being rather small in SN 1997D is an open question. Spectral models appear to rule 
out the case of low metallicity and He enhancement. A possible scenario involves a close 
binary system in which the companion star spiraled in to make the envelope more massive 
( IFodsiadlowski 1992| ). 



The ejection of a small mass of Ni provides a constraint on the mass of the collapsed 
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remnant, M rem ( [Nomoto et al. 1993 ; riiiclcmann et al. 1996 ). In our 26 M Q model, is 
M Tem = 1.8 M Q . If M rom were assumed to be smaller, more 56 Ni would be ejected: e.g. if 
M rem = 1.7 M Q , then M( 56 Ni) = 0.1 M Q . For more massive progenitor models, M rem would 
be larger than 1.8 M & . If the equation of state of nuclear matter is relatively soft, the 
maximum (gravitational) mass of a cold neutron star is ~ 1.5 M Q flBrown fc Bethe 1994|) . 
If this is the case, the remnant of SN 1997D must be a small-mass black hole, and the low 
luminosity tail of SN 1997D might be a signature of the formation of such an object. 
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Fig. 1.— Comparison of the absolute V light curve of SN 1997D with those of SNe 1987A, 
1993J and 1979C. The epochs of maxima have been corrected to J.D. 24500420. Only galactic 
absorption has been considered. The upper limits from the negative observations of Rev. 
Evans (private communication) have been also plotted. The insert shows the bolometric 
light curves of SN 1997D and SN 1987A (Catchpole et al. 1987, 1988). The star marks the 
optical+IR spectrophotometric observation of Jan. 30. 

Fig. 2.— Comparison of the spectra of SN 1997D and SN 1987A. Top: the spectrum of 
SN 1997D at discovery is compared with that of 1987A 1.5 months past explosion. Middle: 
spectra taken 100 days later. The spectrum of SN 1997D is the result of the merging of 
two spectra with different resolutions and signal-to-noise ratios. Bottom: spectra taken 10 
months past maximum. All spectra have been normalized to the peak of the Ha emission 
and corrected for redshift. 

Fig. 3. — The observed absolute V photometry of SN 1997D compared to the model light 
curve (mass = 24 M ; radius = 300i? Q ; mass of 56 Ni= 0.002 M ; energy= 4 x 10 50 ergs). 

Fig. 4. — The observed spectrum of SN 1997D on Jan. 17 (continuous line) and the synthetic 
spectrum computed from the explosion model (dashed line). The parameters of the model 
are: v ph = 970 km s" 1 ; radius= 6000 R Q ; T eff = 6400 K; mass above the photosphere 
M — 14.5 M ; epoch=50 days; solar abundances. 
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